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AN EXPERIMENTAL INVESTIGATION OF THE FREQUENCY 
AND VISCOUS DAMPING OF LIQUIDS 
DURING WEIGHTLESSNESS 

by Jack A. Salzman, Thomas L. Labus, and William J. Masica 

Lewis Research Center 

SUMMARY 

As a part of the general study of liquid behavior in weightlessness, an experimental 
drop-tower investigation was conducted to determine liquid sloshing characteristics in a 
zero Bond number environment. The frequency and damping of the fundamental asym- 
metric sloshing mode in right-circular cylinders (radii ranging from 0. 159 to 1.27 cm) 
were correlated in terms of known system parameters. Contact angles were restricted 
to very near zero degrees such that the sloshing equilibrium interface was highly curved. 
Liquid depths were greater than 1 diameter from the cylinder bottom. 


INTRODUCTION 

A knowledge of the dynamic behavior of liquid propellants in space vehicles has long 
been of fundamental importance to their control and stability. Recently, the reality of 
long-range space flights in a weightless or zero gravity environment has created new 
problems in propellant management. In zero gravity, or, more properly, in a low Bond 
number environment (Bond number indicates the relative effects of gravity to capillary 
forces), positioning of the liquid for venting and restart operations becomes of major 
concern. Investigations establishing the static configuration of the liquid-vapor interface 
in low Bond number environments have been conducted and the results have been doc- 
umented (e. g. , refs. 1 and 2). Actual space mission requirements, however, include 
numerous disturbances which may displace the liquid-vapor interface from its static- 
equilibrium position and induce liquid sloshing. If the frequency of the excitation dis- 
turbance is near the natural frequency of the liquid-vapor interface and especially if the 
damping is small, the slosh amplitude may become large, creating problems in vehicle 



control, liquid level sensing, venting, and restart operations. Hence, the frequency and 
damping associated with the fundamental mode of liquid sloshing in low Bond number 
environments become basic factors in the propellant management problem. 

Investigations of the natural frequency of the liquid-vapor interface in cylinders 
under various system conditions (i. e. , Bond number, contact angle, liquid depth, etc.) 
are presented in references 3 and 4. For the limiting case of a reasonably flat inter- 
face, the theories agree quite well with experimental data. However, as the contact 
angle and Bond number approach zero, such that the interface is highly curved, the 
solutions require numerical approximations. For highly curved interfaces, only limited 
experimental data are available for comparison with the analytic approximations. 

Examples of analytic studies concerned with the viscous damping of oscillatory 
liquid motion may be found in references 5 to 7. The results of experimental 1-g damp- 
ing investigations (refs. 6 to 9) verify the essentials of viscous damping theory except 
that the empirical magnitude of the damping observed was slightly higher than that pre- 
dicted. Damping investigations of a highly curved interface appear to be nonexistent. 

The purpose of this report is to present the results of an experimental investigation 
conducted at the NASA Lewis Research Center on liquid sloshing phenomena in cylindri- 
cal containers under effectively zero Bond number conditions. Contact angles were 
restricted to very near zero degrees such that the liquid -vapor interface was highly 
curved. Liquid depth was greater than one diameter from the cylinder bottom. The 
data obtained on natural frequency and damping characteristics are presented and cor- 
related in terms of known system parameters. Also presented are observations of the 
liquid behavior during slosh including the effects of contact angle variation and large 
initial amplitudes. 


SYMBOLS 


A o 

initial slosh amplitude, cm 

A n 

n 111 cycle slosh amplitude, cm 

a 

2 

system acceleration, cm/sec 

Bo 

p 

Bond number, Bo = aR /|3 

g 

o 

acceleration due to gravity, cm/sec 

k 0 ,k v k' v k 2 

nondimensional constants 

L 

interface height above cylinder bottom 

n 

cycle number 

R 

cylinder radius, cm 



t time, sec 

tj, formation time, sec 

Xj instantaneous interface displacement from sloshing equilibrium at cylinder wall, cm 
a damping coefficient, sec - * 

o o 

j3 specific surface tension, c r/p, cm /sec 
A. damping factor, ot/w Q 

5 logarithmic decrement, 5 = (1/n) In (A n )/ (A n+ jj 

9 contact angle, deg 

fi viscosity, cP 

v kinematic viscosity, j it/p, cm /sec 

3 

p density, g/cm 

o surface tension, dynes/cm 

2 2 3 2 

Cl natural frequency parameter, Cl = (R u>£/ /3) 

damped system frequency, rad/sec 

a> 0 system natural frequency, rad/sec 

u> 0 average measured frequency, rad/sec 

APPARATUS AND PROCEDURE 

The experimental investigation was conducted in a 2. 3-second drop tower (fig. 1). 

-5 

Aerodynamic drag on the experimental package was kept below 10 g by allowing the 
package to fall freely within a protective drag shield. The experiment package which 
was used to conduct the slosh study is shown in figure 2. An air cylinder was used to 
apply a lateral disturbance to a low -friction slide. Precision diameter borosilicate 
glass cylinders, ranging in radii from 0. 159 to 1. 27 centimeters, were mounted on 
this slide. An adjustment of the air-cylinder stroke and pressure to obtain the desired 
initial slosh amplitude was determined experimentally and set prior to each drop. A 
16-millimeter high-speed motion picture camera was used to record time and displace- 
ment characteristics of the liquid behavior. These films were subsequently analyzed on 
a motion analyzer which magnified the image eight times and permitted displacement 
measurements down to 0. 02 centimeter. In an attempt to view the liquid-vapor interface 
more accurately, a backlighting scheme was employed. Basically, this scheme consisted 
of a tapered section of clear plastic with the side facing the camera vapor -blasted. This 
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Figure 2. - Experiment package. 


method ensured uniform illumination by diffusing the entering light. Time measure- 
ments were obtained by viewing a precision digital clock with a calibrated accuracy of 
0. 01 second. 

Analytic reagent grade liquids were used possessing static equilibrium contact 
angles near zero degrees. The properties of these liquids can be found in table I. Con- 
tamination of the glass surfaces and liquids, which could alter the surface tension and 
contact angle of the test liquid, was carefully avoided. The glassware was cleaned 
ultrasonically in a detergent solution, rinsed with distilled water and dried in a warm 
air dryer. 

During the test drop, the application of the lateral disturbance was preceded by a 
predetermined time increment to allow the zero-gravity configuration to form. The 
minimum time required for the interface to reach a sufficiently quiescent configuration 
was empirically determined to be 


j 



( 1 ) 


While this allowed formation time was generally not sufficient to ensure a completely 
quiescent, highly curved interface, the motion was sufficiently damped so that it had no 
effect on the ensuing slosh wave. 
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TABLE I. - SUMMARY OF PARAMETERS 




PREVIOUS INVESTIGATIONS 


It was assumed that, during the fundamental slosh mode, the instantaneous displace- 
ment of the liquid -vapor interface at the cylinder wall can be described by an equa- 
tion of the form 


X t = A 0 e- at sin(w d t) (2) 

The damped system frequency o> d is related to the natural frequency (o Q by the rela- 
tion 


w d = K ' “ } 


.2 „2^ /2 

o 

such that, if the damping coefficient a is small, 

GO j » CtT 


( 3 ) 


( 4 ) 


For convenience, the damping effect is often described in terms of the logarithmic 
decrement 6 where 


6 = 


2 net 





( 5 ) 


The correlation of the natural frequency a> 0 and the logarithmic decrement 6 in 
terms of basic system parameters is the subject of this report. The study was con- 
cerned specifically with small amplitude asymmetric sloshing about a highly curved 
interface (i.e. , zero Bond number and zero degree contact angle) in a cylindrical con- 
tainer. 

Low Bond number sloshing in right-circular cylinders has been investigated in 
references 3 and 4. These linearized inviscid analyses include the effects of contact 
angle and interface curvature on the system sloshing frequency. For the special case of 
a flat equilibrium surface (or a 90° contact angle), the fundamental asymmetric sloshing 
frequency may be expressed nondimensionally as 

n 2 = [6.255 + 1.841 Bo]tanh(l. 841 L/R) (6) 

9 

where the frequency parameter constant S2 is defined as 
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( 7 ) 


0 

For large Bond numbers and liquid-depth ratios, these expressions reduce to the famil- 
iar normal -gravity equation 


J Q = 1.841 g/R (8) 

For other contact angles at arbitrary Bond numbers, no simple relation exists and 
numerical techniques are required for solution. A suggested empirical approximation is 
given in reference 1 of the form 

f2 2 = [6. 255 + 1. 841 Bo - 4. 755 cos 0]tanh 1. 841(L/R) (9) 

For zero Bond numbers, zero degree contact angles and L/R > 2 equation (9) reduces to 

S2 2 = 1.5 (10) 

An extrapolation of the results of Satterlee and Reynolds (ref. 3, in particular, fig. 5) 
and the limited experimental data applicable to zero degree contact angles (ref. 3, p. 178) 
indicates a value for equation (10) approximately between 2 and 3. 

Normal-gravity, high Bond number, investigations of viscous damping in cylinders 
may be found in references 6 to 9. These studies consider small amplitude oscillations 
and implicitly assume 90° contact angles. Their results have shown that viscous dissipa- 
tion in an assumed laminar boundary layer at the cylinder wall and bottom is the primary 
cause of damping. For large liquid-depth ratios (L/R > 2), the damping is practically 
a result only of dissipation at the cylinder wall. Its magnitude is described by the re- 
lation (ref. 8) 


6 


= V 1 / 2 k' 3 / V 1/4 


or its equivalent form 



( 11 ) 


( 12 ) 
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obtained by substituting equation (8) into equation (11). The value of Kj in equation (12) 
for high Bond number sloshing has been empirically determined to be about 6. 1 (ref. 8). 
Investigations of the damping in a low Bond number environment appear nonexistent. 

RESULTS AND DISCUSSION 

Behavior of the Interface During the Fundamental Slosh Mode 

A schematic of the liquid-vapor interface profile during the fundamental slosh mode 
is shown in figure 3. In this study, the height of the lowest point on the interface from 
the cylinder bottom (L) was maintained greater than 2R. The instantaneous displace- 
ment of the edge of the interface at the cylinder wall from the sloshing equilibrium posi- 
tion is defined by . For small slosh amplitudes (Aq ^ 0.5 R), the vertex of the inter- 
face remained stationary with respect to the longitudinal axis of the cylinder. The 
smooth curvature of the interface is representative of the small amplitude, fundamental 
slosh mode observed in all test runs. 

Figure 4 shows a series of photographs illustrating the interface behavior during a 
typical test run. The initial amplitude was produced by first accelerating the experiment 
cylinder laterally (fig. 4(b)) and then providing a quick-stop impulse (fig. 4(c)). Al- 
though the interface is distorted by the large impulse caused by the quick stop (figs. 4(b) 
and (c)), it quickly assumes its well-behaved normal sloshing configuration as can be seen 
in the subsequent oscillations (figs. 4(d) and (e)). 

The motion of along both boundaries was recorded during each test run. A 
representative set of data showing almost two complete cycles of oscillation is given in 
figure 5. In general, for small amplitudes, the displacements of the two edges were 



Figure 3. - Interface profile during fundamental slosh mode. 
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Figure 5. - Interface edge displacement, (luring representative data run. Cylinder radius, 
0.635 centimeter; specific surface tension, 11.8 centimeters cubed per second squared. 


approximately equal and in phase. Although, normally, the position of the interface edge 
is easily measured, at one point in the test run judgement of the actual edge position is 
required. When the cylinder motion is abruptly stopped, the right edge of the interface 
which has been initially displaced upward (fig. 5) continues to advance while the re- 
mainder of the interface begins to shift toward the opposite boundary. The initial ad- 
vancing edge then thins out until only a film of liquid remains. The actual interface edge 
gradually reappears below the equilibrium position. Thus, there is no sudden drop in 
the edge position as might be assumed from the initial portion of figure 5, but rather 
only a lapse in its determination. The residual liquid film left on the cylinder wall 
caused by the initial amplitude produces a slight drop in the sloshing equilibrium position 
as compared with the zero Bond number equilibrium configuration (e.g., fig. 5). No 
additional change in the sloshing equilibrium position was measured during the remaining 
test time. The liquid film remains on the wall and sloshing actually occurs on a thin 
layer of liquid rather than on an unwetted surface. 


Natural Frequency Results 

The half -periods of oscillation in each test run were measured and the corresponding 
frequency values were calculated using the recorded positions of both interface edges. 

A tabulation of these frequencies and the subsequent average frequency of each run is 
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Figure 6. - Experimental correlation of natural frequency. Spe- 
cific surface tensions, 11. 8 to 29. 9 centimeters cubed per 
second squared; cylinder radii, 0. 159 to 1. 270 centimeters. 


listed in table I. The periods of oscillation of any particular system were relatively 
constant as indicated by the frequencies listed. From these data, the empirical value 
for Q, in equation (7) was determined to be 2. 6 by calculating its value for each mea- 
sured frequency (i. e. , each half period) and then computing a simple arithmetic average. 
As a measure of uncertainty, the standard deviation of was calculated and deter- 
mined to be 0. 41. A plot of the average frequencies as a function of the ratio 
is shown in figure 6 along with the calculated line representing the equation 

Wq = 2. 6 o/pR 3 (13) 

From this figure it can be seen that the experimental data does follow the relation given 
by equation (13). 

It is realized that the measured frequency shown is not truly the natural frequency 
of the system, but rather a damped frequency. However, the damping coefficient was 
sufficiently small such that the approximation in equation (4) was valid. Calculations 
using the measured values of 6 showed that the variation was less than 5 percent. 

From equation (1) which predicts formation time requirements, and equation (13), 
which yields the time required for a period of slosh oscillation, it is evident that the 

Q 

limitations on the lower value of j3/R and consequently the highest value of R used in 
this study were imposed by the time restrictions of the test facility (i. e. , 2.22 seconds 
of zero-gravity time). Values of /3 ranged from 11. 8 to 29. 9 cubic centimeters per 
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second squared; R values were varied from 0. 159 to 1.270 centimeters. 

Contact-angle effects . - All liquids employed in this study exhibited near zero 
degree static-contact angles on the cylinder surfaces. Since the contact angle will affect 
the natural frequency (refs. 3 and 4), particular care was taken to avoid inducing 
finite -contact angles during the slosh motion. As indicated in reference 10, finite - 
contact angles could be introduced during flow, even though the static angle is zero. 
While the results of reference 10 are strictly applicable only for steady flow, the effect 
of. viscosity as reported was of equal importance in sloshing motion - the larger the vis- 
cosity, the more likely a finite -contact angle. The photograph in figure 7 shows the 
large contact angle induced for a fluid with a viscosity of 15 centipoise. No measurable 
contact-angle variations were observed for the data in table I. However, when using 

butanol (viscosity of 2. 9 cP), it was easy to obtain large contact-angle variations. The 

o 

butanol data were not used in determining because of the possibility of erroneously 
raising its value. 


Viscous Damping Results 


Damping of the interface oscillations was determined by observing the decay of the 
slosh amplitudes A n of each successive cycle. As an example, the amplitudes re- 
corded in figure 5 are plotted in figure 8. On this semilogarithmic plot, the straight 
line representation of the data indicates that the damping is exponential. This was 
representative of all small amplitude test runs. Therefore, the damping magnitude may 
be measured in terms of the logarithmic decrement defined in equation (5). The magni- 
tude of 6 was obtained by hand -fitting a straight line through the points (e. g. , fig. 8). 
The values of 6 are tabulated for each data run in table I. 

Dimensional analysis of the logarithmic decrement and the pertinent system vari- 
ables yielded the grouping, p/poR. This grouping together with a parametric analysis 
of the data yielded the following damping relation 



(14) 


The value of K 2 was calculated from each data point, and its average value was deter- 
mined to be 21.6 with a standard deviation of 2. 2. A plot of equation (14) along with the 
data is shown in figure 9. Equation (14) may be rewritten in terms of the natural fre- 
quency permitting comparison with normal gravity damping by substitution of equa- 
tion (7) which yields 
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Logarithmic decrement, 6 




Figure 10. - Zero Bond number damping as a function of system variables and meas- 
ured frequency. 
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1/2 


(15) 


6 =K 



The forms of equations (12) and (15) are similar; this would tend to indicate that the 

cause of damping in zero gravity is viscous dissipation in a boundary layer adjacent to 

the cylinder wall. A plot of equation (15) along with the data appears in figure 10. 

Using the measured values of u> 0 and 6 yielded a value of K'^ of 28. 1 with a standard 

deviation of 2.5. This value also agrees with the value of K', obtained when the in- 

2 1 

dependently determined constants £2 and Kg are used. Comparing K’j and Kj in 
equation (12) for high Bond numbers indicates that the constant is about four times 
greater in a zero Bond number environment. 


Amplitude Dependence 

The empirical correlations describing the natural frequency and logarithmic 
decrement are valid for slosh amplitudes less than 0.5 R. In this amplitude range, the 
natural frequency was observed to be constant for any particular system (i. e. , it did 
not vary with amplitude) and the damping was observed to be exponential. An attempt 
was made to determine the maximum amplitude to which the frequency and damping cor- 
relations remained valid, but no conclusive results were obtained. While it was 
possible to obtain initial amplitudes as high as 1. 0 R, the amplitude after the first 
half cycle of oscillation always decayed to less than 0. 5 R. Since several half periods of 
oscillation were required for a frequency measurement, the effect of amplitudes in this 
range on frequency could not be obtained. It is noted, however, that, for amplitudes 
above 0. 5 R, the initial amplitudes were generally not equal and the position of the 
vertex of the interface was no longer stationary. 


SUMMARY OF RESULTS 

An experimental investigation was conducted to determine the natural frequency and 
damping characteristics of the liquid-vapor interface in a zero Bond number (zero- 
gravity) environment. The study employed flat-bottomed right-circular cylinders and 
and liquids possessing near zero degree contact angles on the cylinder walls. The 
liquids had viscosities /i between 0. 32 and 2. 9 centipoise and specific surface tensions 
/3 between 11.8 and 30. 4 cubic centimeters per second squared. The range of radii R 
employed was from 0. 159 to 1. 270 centimeters. Under the stipulation that the liquid- 
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vapor interface is greater than 2 radii from the tank bottom, that the initial slosh 
amplitude is less than 0. 5 R, and that the contact angle does not vary during sloshing, 
the investigation yielded the following results: 

1. The constant O in the natural frequency relation, w Q = £1 /3/R , was empir- 
ically determined to be 2. 6, where co Q is the natural frequency. 

2. The damping in terms of the logarithmic decrement (6) was experimentally 
determined to follow the relation 6 = K 2 (m 2 /p°'R)* // 4 where Kg was empirically deter- 
mined to be 21.6 and p and a are liquid density and surface tension, respectively. 

3. The logarithmic decrement was also expressed as 6 = K^ju/pu^R 2 ) 1 / 2 permit- 
ting comparison with normal gravity results. K’j was determined to be 28. 1, which is 
approximately four times greater than Kp the normal gravity constant. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 4, 1967, 

124-09-03-01-22. 
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